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I. INTRODUCTION 



We consider the dynamical behavior of a quantum system governed by a time dependent 
Hamiltonian H{t), characterized by the following properties: 

(a) H{t) is a smooth family of self-adjoint, bounded operators, and 

(b) the time derivative H{t) is compactly supported in the interval [0,r]. 

One can think of such a family as a switching system, i.e., a system that coincides with 
H{t) = Hi in the past (t < 0), and switches to the system H{t) = Hp in the future (t > t). 

Our goal is to establish an upper bound on the minimal running time r needed to make 
the error small in the adiabatic theorem for the switching Hamiltonian H[s). Our results 
substantiate the ideas presented in Ref. [sl. 

It is convenient to rescale the time t to s = t/r. With a minor abuse of notation, the 
system then evolves according to the Schrodinger equation: 

iMs) = tH{s)Ms), MO) = i^i, (i-i) 

where ipi is the ground state of Hj. The adiabatic theorem of quantum mechanics ensures 
that under certain conditions, if ipri^) is close to the ground state of Hj, then ipr{^) is close 
to the ground state of Hp (see the theorems below for details). 

Recently it has been realized that the adiabatic approximation could be used as the 
fundamental ingredient for a method of quantum computation^. This adiabatic quantum 
computing (AQC) has generated a resurgence of interest in the adiabatic theorem. In gen- 
eral, quantum computing attempts to exploit quantum mechanics to obtain a speedup for 
classically difficult computational problems. It was subsequently shown in Ref. [ij that AQC 
provides a universal model, equivalent (in terms of complexity) to the quantum gate model 
(and hence to other universal models). In adiabatic quantum computing, one solves a com- 
putational problem by using an adiabatically changing Hamiltonian function whose initial 
ground state encodes the input and whose final ground state encodes the output. The time 
T taken to reach the final ground state is the "running time" of the quantum adiabatic algo- 
rithm. One would like to minimize it, while at the same time keeping the distance between 
the actual final state and the desired final ground state small. The crucial parameter on 
which r depends is the minimal value g of the spectral gap g{s) between the ground state 
energy of H{s) and the rest of its spectrum. 



Adiabatic theorems fall into two categories: those that describe the solutions for all 
times, including times s G [0, 1], and those that characterize the solutions only for large 
times s > 1, where the Hamiltonian is time-independent again. Interestingly, the latter give 
more precision for long times. We call the first category, the one that applies to all times, 
"uniform" ; the second is the "long time" category. 

A representative result from the uniform category is the following: See, e.g., Ref . [2I: 

Theorem (Uniform adiabatic theorem). Suppose H{s) is a t -independent, twice differen- 
tiable family of bounded self-adjomt operators on the interval [0, 1]. Suppose m addition 
that 

g := min g{s) > for all s G [0, 1] . (1.2) 

se[o, 1] 

Then, for any s G [0, 1], the solution tpris) to the initial value problem f ll.l|) satisfies 

dist {^r{s), Range P{s)) = O (r^^) , (1.3) 

where P{s) is the orthogonal projection onto the corresponding eigenstate of H{s). 

A characteristic result from the long time category is 

Theorem (Long time adiabatic theorem). Suppose H{s) is a t -independent, C°° family of 
bounded self-adjoint operators that satisfies (II. 2|) . If H{s) is supported on [0, 1], then the 
solution iPt{s) to the initial value problem (11. 1|) satisfies 

dist (ipris), Range Pp) = o for all s>l, (1.4) 

for any n E N. 
Remarks: 

1. One can summarize these results by saying that slowly starting and finishing the 
interpolation decreases the error. 



2. In general, there is no uniformity in n in ( 11.4^ : the term on the right hand side is of 
order c„r~" where c„ grows rapidly with n {c.f., the following discussion). 

3. The distinction between the uniform and long time adiabatic theorems has an analog 
for integrals. Suppose g{s) G C°°(M) has support in [0, 1]. Then 



g{t) e'^^dt 







o(r-") if s > 1 ; 
O(r-i) if s G (0, 1) . 



3 



To describe our result we begin by introducing some notation. Let Hj and Hp be two 
self-adjoint operators on a Hilbert space H that satisfy ||-f^/|| = = 1- Let H{s) be 

a C°° family that switches between Hj and Hp as above, and let P{s) be the orthogonal 
projection onto the eigenvalue E{s) of H{s). 
We assume the following hypothesis: 

Assumption I.l (Minimal gap). For all s G [0, 1], we assume the operator H{s) has 

an eigenprojector P{s), with eigenenergy E{s) separated by a gap g{s) from the rest of its 

spectrum. We assume g = min g(s) is strictly greater than zero. 

se[o, 1] 

Remarks: 

1. The eigenprojection P{s) is allowed to be degenerate. In particular, it can be infinitely 
degenerate. 

2. If Rank P/ < oo, then it follows from analytic perturbation theory that 
Rank P{s) = Rank Pi for all s e [0, 1]. 

In the quantum computational setting the typical value of the gap g is very small. We are 
interested in minimizing the running time r so that the error in the adiabatic approximation 
is small for such values of g. We therefore investigate how the coefficients c„ depend on the 
gap g, and we minimize the running time r so that CnT^" = o(l) for an optimally chosen 



20, which states that the running time 



value of n. One recent result in this direction is Ref. 
T of the order g~^ makes the error in the long time adiabatic theorem small. 

It has to be noted that as long as one is interested in just making the error in the 
adiabatic theorem o(l) in (7 <^ 1 rather than say 0{g^), there is not much difference between 
the uniform and long time adiabatic theorems. Our main assertion below is consequently 
formulated as a uniform adiabatic theorem: 

Theorem 1.2. Under Assumption \I.l\ the error dist (ipris), Range P{s)) for s G [0, 1] in 
the adiabatic theorem is o(l) for g ^ I, whenever 



T > Kg ^ I In (7 



1 6« 



for some g -independent constant K > 0, if the Hamiltonian belongs to the class G" with 
a > 1, given in the following definition. 
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Definition 1.3. An operator valued function H{s) belongs to the Gevrey class C^lR), Ref.\^, 
if H{s) is supported in the interval [0, 1] and there exists a constant C, such that for any 
k>l, 



max 

sG[0,l] 



df'i 



We define = |J 



R>0 



When a = 1, this class coincides with the set of analytic functions, and the only such 
functions are constants. For a > 1 there are functions in the class that are not constant. 



A. A Prototypical Example: An Interpolating Hamiltonian 

We call H{s) an interpolating Hamiltonian if 

H{s) := il-fis))Hj + fis)HF, 
where / is a monotone increasing function on M that satisfies 

/eC°°(M); supp/c[0,l]; /(0) = 0; and /(I) = 1 
Specifically, we can construct / G C°°(R) as follows: 



(1.5) 



(1.6) 



f{t) = / g{s)ds, where g{s) 

J —oo 







if s^[0, 1]; 



/3 exp (- ^) if s G (0, 1) . 



Here, /3 is a normalization constant, chosen so that /(I) = 1. For this family, we have 

d'^s d'^s 

Hence, H{s) G G^. 
Remarks: 

1. For analytic families of Hamiltonians one can obtain sharper control of the transition 
probability in the adiabatic approximation. For two level systems this goes back to 



Landau and Zener, w 
Refs. 
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lo showed that the transition probability was 0(e '"^^'^); see 



14, and 



16l for rigorous treatments. The analogous statement (albeit 



with the less explicit dependence on ^f) is known to hold for the general analytic 
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H{s); see Ref. |22| and references therein. For the robust adiabatic theorem with 
Rank i^i? ^ dimTi, one can develop a lower bound on the run time of the form 
r = 0{g-y\\ng\), see Ref.L 



2. The above remark shows that our result is nearly optimal. In the switching family 
setting, the requirement that H{s) have compact support exacts a price on the shortest 
achievable run time r. In our result, it introduces logarithmic corrections to the 
Landau-Zener type bound r = o^g'"^). 

3. When the gap g{s) becomes small only for finitely many times s, such as in the Grover 
search problem^, one can devise a gap-sensitive interpolating function f{s) that yields 



a much better estimate, t = 0{g ^), see e.g., Ref. 



13 



and references therein. 



B. Relation with Past Work 

Mathematical analysis of adiabatic behavior has a very rich history starting with the first 
rigorous result by Katc^* for rank one projections P{s) and Nenciu22, for more general P{s). 
We do not attempt to give an exhaustive survey of the related literature, but rather focus 
on articles that steadily improved the understanding of the long time adiabatic theorem. 

The first mathematically rigorous result in this direction goes back to Lenardr^. He proved 
a long time theorem for a Hamiltonian function of finite rank with no level crossings. The 
next significant progress came with the 1987 work of Avron, Seller, and Yaffe^. They proved 
a long time result for a general family of Hamiltonians. They did not consider dependence 
on the gap g. In 1991 Joye and Pfister— obtained an estimate on the exponential decay rate 
for the 2x2 matrix case. Three years later Martinea^ realized that the adiabatic transition 
probability could be considered as a tunneling effect in energy space. He used microlocal 
analysis to prove an analogous estimate for the general case. In 2002, Hagedorn and Joye^i 
proved exponential error estimates for a special case using only elementary techniques. In 
2004, Nakamura and Sordoni^ combined these tunneling estimates in energy space with the 
stationary theory of time-dependent scattering to simplify the method of Martinez. Both 
the Martinez and Nakamura-Sordoni results depended on the choice of a point E{s) in the 
spectral gap for H{s). Martinez assumed E{s) was analytic in s. Nakamura and Sordoni 
assumed it was in the Gevrey class and satisfied lE^'^^s)] < (1 + s'^)~p for some constants 
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Ck and p. Gevrey class Hamiltonians were also considered in the unpublished works of 
Jungi^ and Nenciu^^. 

In this work we propose a straightforward method for computing an upper bound on the 



running time, based on Nenciu's expansion, introduced in Ref. l2J. Our method might yield 
less precise estimates than those obtained by microlocal analysis. Our estimate contains the 
logarithmic prefactor, but our results do not rely on analyticity of E{s). 



We are grate 
purpose of Ref. 



ul to the referee for making us aware of the unpublished preprint^. The 



25 



was to quantify the (/-dependence of the error estimates for the long 
time adiabatic theorem in Ref. I2J in the g ^ 1 limit. Although the asymptotics for (7 ^ 1 
established in the present article are different from those of Ref. 



25 



Ref. 



, the basic strategy of 



25l is very similar to ours. 



C. Proof strategy 

The analysis of the wave function iPt{s) is hampered by the fact that it carries the highly 
oscillating, memory-dependent phase. In particular, one cannot easily decompose it into an 
asymptotic series in powers of r^^. However, the projector Pt-{s) onto the state ipris) (given 
by Pt{s) = \ipr{s)) {'4't{s)\) naturally has no phase, and admits the asymptotic expansion of 
the form 

Pr{s) ~ Bo{s) + -B,{s) + \b2{s) + .... (1.7) 

See the next section for details. 

A curious fact about this expansion is that it is instantaneous in H{s) - no memory 
term is present in any finite order r^K The first term in (11.711 satisfies Bq{s) = P{s). 
Moreover, Bj{s) = for all j > 1 provided that all derivatives of H{s) vanish. In particular, 
Pt{s) — P{s) = o(r^") for any n E N and s > 1, at least formally. To make the argument 
rigorous, one has to estimate the remainder Rj^ of the series (after the truncation at some 
Bn). It turns out that H-RatH is proportional to r^^ and depends linearly on B^lt), for 
t G [0, s] (the remainder encodes the memory of the process). Therefore one wants to 
investigate the dependence B^it) on the size of the gap g{t) and on the integer N. We 
shall see in Section UTIl that in general, ||-BAr(t)|| < CNg~^'^'^^^\t), where Cn grows (super) 
factorially fast in A^. The proof of the uniform adiabatic theorem. Theorem [H usually 



involves truncation of the expansion (JLTj) at A^ = 1, with the consequent estimate on the 



error of the form Ri = 0{T~^g~^). For a smooth family H{s), however, one can use the 
fact that the factor r^'^ g~^'^^+'^) decreases with N whenever r^^g^"^ < 1. Therefore one can 
look for the optimal value A^^opt for the truncation of the asymptotic series. Evaluation of 
||i?7Voptll ^^'i estimation of the sum of the first A^opt terms in the expansion yields Theorem 
11.21 (see Section HV] for details). The logarithmic correction in the theorem is an artifact of 
the appearance of the extra factor g~^ in the expression g~^'^^^^\ 



II. AN ASYMPTOTIC EXPANSION FOR 



We now consider the initial value problem for the Heisenberg equation, 

= r [His) , Pris)] 
= Pi 




(II.l) 



In 1993, G. Nenciu^i found a general form for the solution to this Heisenberg equation. 
His idea was to look for an asymptotic series of the form (II. 7p . Substituting (II. 7p into pi.ip . 
we obtain a sequence of differential equations 



iiS^is) = [H{s),Bj+,{s)], J = 0,1, 



(II.2a) 



In addition, since -Pt(s) is a projection for each s, we have Pt{s)'^ = Pt{s)- This generates 
the following sequence of algebraic relations: 



Bj{s) 



E 

m=0 



Bm{s) Bj_rn{s), j = 0, 1, 



(II.2b) 



In particular: i?o(s)^ = -Bo(s), so Bo{s) is a projection for each s. 

It turns out that the system of hierarchical relations ( III.2ap and (III.2bp has a unique 
solution, which is given by the following recursive construction: 



Bois) = P{s) 
where Rz{s) = {H{s) — z)~^, 



P{s),B,^,{s) 



R,is)dz + Sj{s) - 2P{s)Sj{s)P{s) 



(11.3) 



(11.4) 



m=l 



and the contour F encircles only the ground state energy. In particular the first order term 
is given by 



P{s) , Pis) 



RJs) dz . 



(11.5) 



One can truncate the expansion^ 01. 7p at some finite order A; > by observing that 



Pris) = Bois) + -Biis) + 



1 

7n 



1 



r 



N 



Ur{s,r)BN{r)Ur{r,s)dr. (11.6) 



where UT-{s,t) is the unitary Schrodinger propagator that satisfies 



^i-sUris, r) 



Urir, r) 



TH{s)Ur{s,r) 
1. 



(11.7) 



III. ESTIMATES ON 5„ 



Without loss of generality we assume that the constants C, R in Definition 11.31 and the 
value of l/g{s) are all greater than or equal to 1. To estimate the minimal run time r, we 
use the following result: 



Lemma III.l. Suppose H{s) belongs to the Gevrey class C^^R). Then, 



d'B^js) 
dJ's 



< 



1 



[1072 + 0.3)2 



g-2n-k (2Ci?(A: + 3n)2°) 



2a \ fc+3n 



(III.l) 



Proof. We use induction to prove this lemma. 
We define 

1 



L(n, k) 



g-2n-k (^2C Rik + Sn)^") 



k+3n 



(10 n + 0.3)2 

Since Bq{s) = P{s), we have ||i?o(s)|| = 1. To estimate Bo{s), we use the representation 



Bo{s) 



Rz{s) dz , 



where Rz{s) := {H{s) — z)~^ and the contour F is given by 

{zeT: \z-E,{s)\=g{s)/2} 

We note that the circumference of F is it g{s). 
Differentiating both sides of flIII.2p and using 



(III.2) 



Rz{s) = - Rz{s)H{s)Rz{s), 



we see that 



27ri 

Since 



Bo{s) = - ^ / Rz{s)H{s)R,{s)dz. (III.3) 



r 



and our assumptions require 



\Rzis)\\ = ^ for zer (III.4) 



\H{s)\\ < CR, (III.5) 



we obtain 

\\Bo{s)\\ < C^^. (III.6) 

Recall now the definition of an integer composition: If /c is a positive integer, then a 
composition vr^ ; of k is an ordered set of positive integers Pi, P2, ■ ■ ■ , Pi whose sum is k. Let 
Uk denote the set of all possible integer compositions Hk,i of k, and let 11^;; C 11^ denote the 
set of compositions of k into exactly I parts. We note that the cardinality of 11^ is 2''~^. 

Armed with this notion, we observe that 

d'^RJs) v^, ,w V- TT fr. . N dP^H(s)] \ ^ , , ,,,, 

1=1 TTk^ieHk^i yPiETrfc,, ^ y 

where the Cj^^ ; are the multinomial coefficients 



Pu---,Pi 



In particular. 



by the multinomial theorem. 

Taking the norm on both sides of 0111.7^ and using equations IIII.4I and IIII.9I as well as 
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the fact that a > 1, we obtain the following bound: 



dJ^s 



1=1 



l+l 



(P^H{s) 
dP^s 



^ E c" (-^) E n 

4 f 2CRk'^'^V 



This initiates the induction since 



(s)|| < g/2 sup <2 



2CRk 



2a \ k 



9is) 



< L{0,k) 



Next, we verify the induction step: Suppose 

d'^BJs) 



< L{n, k), 



for some n and all k. 

For z G r, we use definition ( 111.3^ and the Leibniz rule to bound 

d'Bn+i{s) 



d^s 

<9{s) E 



k\ 



ki+k2+ 



ki\k2\kz\ki\ 



+ 2 E 



k\ 



k-i+k2+k-j,=k 

d'^s 



d^\s 

d'^'Bois 
d^^s 



d^^s 

d^'Sn+l{s 

d^^s 



d^'^^'Br.is) 
d^^s 



(III.IO) 



d^^RA 
d'^^s 



(III.ll) 
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Using f lIL4p and the induction estimate (lIII.ip . we bound 



< 



E E 



hi hi 



d'^^s 



d'^^Bn+l-ijs 

d^^s 



i=l ki+k2=k 

_ k\ ih + 3^)'" + 3(n + 1 - (te+3(n+i-i)) 

k 

i=l k\+k2=k 

0.05 



< 9' 

n 

E E 



< 



(10n + 10.3)2 
0.05-L(n + l, k) 



ki\k2\ (10i + 0.3)2 (lOn + 10.3- 10z)2 

g-2in+i)-k (2Ci?)'=+'("+^) (A; + 3(n+l))'"('=+'("+^» 



where, in the last inequahty, we have used Corollary IV. 31 
We use this bound together with (1III.10|) to verify that 



< 



2 E 

0.4 

(10n+ 10.3)2 



k\ 



fci! k2\ ksl 



d^^s 



d''Sr,+l{s) 



d'^^s 



g-2in+l)-k (2Ci?) 



fc+3(n+l) 



E 



k\ 



d'^^B, 



d^^s 



k-i+k2+kz=k 



k,\ k2\ k^\ 



(A;i)2°^-i {k2 + 3(n + i))2"(^'2+3(n+i)) (^3)2^/03 



< 



0.4 



(10n+ 10.3)2 
0.4-L(n + 1, k) 



g-2(n+l)-k (2Ci?)'=+3('^+^) (fc + 3(^ + l))2"(fc+3(n+l)) 



(III.12) 



where, in the last inequality, we have used Lemma IV. 41 
Next, we bound the first contribution in fllll.lip . 
We first observe that for n = k = 0, we have 



9 \ms)\\ \\Bois)\\ 



dBojs) 
ds 



\\RM\\=4g-' {2CR) 
< 0.02 ■L(l, 0). 



(III.13) 
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For n = and k = 1, we have 



E 

fcl+fc2 + 
^3+^4=1 



d^''' + ^Bn{s) 

d^--^+^s 



< T' 2'" g-^ {2CRy 

< 0.001 ■ 1) . 



(III.14) 



For 2n + k > 2, we bound 



ki+k2 + 
k3+k4=k 



k\ 



ki\ k2\ ksl k^\ 



d'^^s 



d^'Bois) 
d^-^s 



d'^-'+'B^s) 



d^^s 



< 



(lOn + 0.3) 



fc+l+3n 



X 



ki+k2 + 
k3+k4=k 



k\ 



{k^T'^' (fc3 + l + 3n) 



A;i! fca! /ta! ^4! 



2a(fc3+l+3n) ^^(^^^"^4 



< 2 



0.3 



(10n + 0.3)2 
< 0.5 -1.(71+1, fc). 



g-2{n+l)-k (2Ci?)'=+'+='" (fc + 1 + 372) 



20 {fc+l+3n) 



(III.15) 



where we have again used Lemma IV. 4[ 
Combining ( IIII.lip - OIII.ISP we arrive at 



d^Bn+i{s) 
d^s 



< L(77+ 1, k). 



This proves the lemma. 



□ 



IV. PUTTING EVERYTHING TOGETHER 



To estimate the error ||-Pr('S) — -P(s)||, we use (III.6p . which yields 



^ 1 

Pr{s)-P{s)\\ < 11^. 



-N 



— / \\BN{r) II dr 



(IV. 1) 
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To find the optimal value for A^, we estimate the last term first. Using the bound f llll.ip . 
we see that 



1 



5jv(r) II dr < (2 C + SiV)^-) 



(r/g) 



1/3 



2(:7i?(3A^ + 1)2° 



(IV.2) 



(^^2)1/3 

We want to find the value of that minimizes the right hand side. Differentiating, we see 
that the minimizing value of satisfies 



^2)1/3 y/2« 



2CR 



> 3iVopt + l > 



2CR 



Substituting Nopt into flIV.2p . we get the following bound 



I < {rlgfl^ exp (^-2e-a ( j 

For ^ 1, the above expression is o(l) for small g provided r satisfies 



r > Kg-'\\^gf'^, 



(IV.3) 



for some sufficiently large, ^^-independent constant K > 2. 

To estimate the sum of the first = A'opt terms on the right hand side of (1IV.ll) we again 
use the bound (lIII.ip to get 



Nopt Nopt 
Nopt 



1 



Vrg 



3j 



N 



opt 



1 



17- r 



N. 



1/2 

opt 



for r that satisfies ( 1IV.3P and g sufficiently small. This proves our main result. Theorem II. 2[ 

□ 
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V. APPENDIX 

In this appendix, we collect several technical results. 

Lemma V.l. For a// a > 1, fci = 0, 1, 2, ■ ■ ■ , A;2 = 0, 1, 2, ■ ■ ■ , and = 1, 2, 3, ■ ■ ■ , 

u 

(A;i A;^'^^ < 4-("~^)'^''^('^i+"''=2) (A; + n)" , (V.l) 



(A;i!)(A;2!) 
where k = ki + 

Proof. We first prove the above inequality for a = 1. We begin by noting that 

i-l k^ - - 

{ki + n)''^ k^^ < ^ \ {ki + n)''^ k^' = {k + n)^ . 



Since {ki + n)" < {k + n)", the inequality is preserved if we multiply the left hand side by 
{ki + n)" and right hand side by {k + n)". Hence 



max 



{ki + nf'^'' k^' < {k + n)'=+" . (V.2) 



ki+k2=k (kil) {k2\] 
This proves the result for a = 1. 

We now consider the case a > 1. For < a; < 1, we have x{l — x) < 1/4. Applying this 

. , ki + n , 

with X = , we see that 

k + n 

(ki + n) k2 < ^ {k + nf. (V.3) 

Suppose k2 < ki + n. From (|Vj), we see that {ki + nf'' k\^ < A'^^ (k + nf^\ We 
multiply the left hand side of this inequality by {ki + 72)^=1 +"-^=2^ ^t^^ multiply the right 
hand side by the larger or equal quantity {k + n)^^^""'^^. We conclude that 

(fci + k^^ < A-'"' {k + (V.4) 

Similarly, if ki + n < ^2, then (|VJ]) implies (fci + < 4-('=i+") {k + 

We multiply the left hand side of this inequality by k^^^^^^" and the right hand side by the 
larger or equal quantity {k + 7j,)^2-fci-n_ conclude that 

{ki + nf^^'' k^^ < 4-('=i+") {k + n)^+". (V.5) 
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Combining f lV.4p and (IV. SP and raising the result to the {a — 1) power, we see that 



This and f lV.2p imply the lemma. □ 

Corollary V.2. Suppose a, ki, k2, k, andn are as in Lemma \V.l\ For 1 < i < n, we have 

— {ki + n + l- {k^ + 



h\k2\ 

< 4-(a-l)mm(A:i+l,fc2+l) _|_ ^ _|_ -|^^a(fc+n+l)^ (V.6) 

Proof. Suppose first that ki+n + l — i > k2 + i. Then we have ki + n>k2 + 2i — l>k2 + l- 
For a > 6 > 0, the function g{j) := (a+j)"*^"^-'^ (6 — j)"^''"-'^ is non-decreasing on the interval 
^ j ^ b. We apply this with a = ki + n + 1 — i, b = k2 + i, and j = i — 1. The inequality 
9(0) < gU) yields 

Since k2 + I < ki + n, Lemma IV. II with A;2 + 1 in place of k2 implies 



fcl!(fc2 + l)! 

Therefore, 

(A; + l)! 



fci! (A;2 + 1)! 

< 4-("-i){fe+i) (A; + n + 1)" (^+"+^) . (V.7) 

Suppose next that ki + n+l—i< k2 + i. We again use ^'(0) < g{j), but with a = k2 + i, 
b = ki + n + l— i and j = n — i. This yields 

We note that ki + n + 1 — i<k2 + i implies ki + 1 < k2 — n + 2i. Since i < n, we have 
fci + 1 < /c2 + ^- We apply Lemma IV. II with ki replaced by ^2 and ^2 replaced by ki + 1. 
This yields 
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Therefore, 

(k + iy. 



ih + 1)! h^. 

< 4-("-i)('^i+i) {k + n+ 1)" (^+"+1) . (V.8) 

Inequalities (IV. 711 and ( jV.811 imply ( 1V.6I) because 

A;! . f + (fc + 1)! 

< mm ' 



h\ k2\ - I ki\ {k2 + l)V {ki + l)\ k2\ 



We also have the following consequence of Corollary IV.2I 
Corollary V.3. For a > 1, k = 1, 2, 3, ■ ■ ■ , and n = 0, 1, 2, ■ ■ ■ , we have 



□ 



(10 72+ 10.3- 10 i)2 (10z + 0.3)2 

1=1 ki+k2=k ^ / \ / 

kil k2l 

^ 0.05 ■ 4-(""3/2) (A; + 3(n + l))"(^+3("+i)) 

- 1 _ 4-(a-i) (10n + 10.3)2 ■ ^ ■ 



Proof. Using Corollary IV.2t we estimate the sum over ki + k2 = k hy 

k\ 
kjk^. 



(A;, + 3(n + l-z)r('=i+3(n+i-0) (^k2 + 3tr^'''+''^ 

k\+k2=k 



< (A; + 3 (n + 1))°^''+^^"+-^^^ ^ ^-(a-l) mm(A:i+l,A:2+l) 

ki+k2=k 



< (A; + 3(n + l))°(^'+3("+i)) 2 ■ 4-("-i) ^ 4 

fci=0 

4-(a-3/2) 

= (fc + 3(n + l))"('=+=^("+i)) 



oo 



1 _ 4-(l-a) ■ 

The result follows from this and the inequality 

^ (lOn + 10.3- 10i)2 (10i + 0.3)2 ^ (10 n + 10.3)2 ' *^ ' 

For n = 1, 2, 3, we check this inequality by explicit computation. For n > 4, we first 
note that 

+ I • (v.ii) 



{a — b) b a \ a — b b 
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We take the square of both sides and use a = 10 n + 10.6 and b = 10 i + 0.3. We then sum 
over i from 1 to n. Using the i ^ n + 1 — z symmetry, this yields 

n 

E 



(10n+ 10.3- 10i)2 (10i + 0.3)2 

'A 1 A 1 

(10n + 10.6)2 (10i + 0.3)2 + (10n+ 10.3- 10^(10^ + 0.3) 

El 4 1 



(10n + 10.6)2 j-^ (10i + 0.3)2 (10n + 10.6)3 j-^ (lOz + 0.3)' 

where in the second step we have used (IV. lip and the i ^ n + 1 — i symmetry one more 
time. 

The first term on the right is bounded by 

2/100 ^ 1 1/100 7r2 



(10n+ 10.6)2 ^ ^2 (10n + 10.6)2 3 
The second term is bounded by 

4/100 1 >A 1 

(10 n + 10.6)2 n+l ^ ~i' 



For n > 4, 



— t-<-- 

n + 1 ^ 2 ~ 3 

i=l 



Combining these results, 



1 ^ (7r2 + 5)/300 

^ (lOn + 10.3- 10z)2 (10i + 0.3)2 - (iOn + 10.6)2 

0.05 

< 



;iOn+ 10.6)2 ' 

which proves (IV.lOp . □ 
Arguments similar to the proof of Lemma IV. II generalize to prove the following: 
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Lemma V.4. For a> 1, define k, 



4-(c«-3/2) 



. Then for A; = 1, 2, 3, ■ ■ ■ 



and 



« = 0, 1, 2, ■ . . , 



(V.12) 



fci! k2\ k^\ 



^ k\k\k\k\ + nr^'^^"^^ {k^r'^ {k,r'^ {k,r'- < 4 (k+nr^'+^^K (v.is) 
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